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ABSTRACT 

One of the key predictions of the merger hypothesis for the origin of early-type (elliptical and lenticular) 
galaxies is that tidally-induced asymmetric structure should correlate with signatures of a relatively young 
stellar population. Such a signature was found by Schweizer and Seitzer (1992; AJ, 104, 1039) at roughly 
4cr confidence. In this paper, we revisit this issue with a nearly ten-fold larger sample of 0.01 < z < 0.03 
galaxies selected from the Two Micron All-Sky Survey and the Sloan Digital Sky Survey. We parameterize tidal 
structure using a repeatable algorithmic measure of asymmetry, and correlate this with color offset from the 
early-type galaxy color-magnitude relation. We recover the color offset-asymmetry correlation; furthermore, 
we demonstrate observationally for the first time that this effect is driven by a highly- significant trend towards 
younger ages at higher asymmetry values. We present a simple model for the evolution of early-type galaxies 
through gas-rich major and minor mergers that reproduces their observed build-up from z = 1 to the present day 
and the distribution of present-day colors and ages. We show using this model that if both stellar populations 
and asymmetry were ideal 'clocks' measuring the time since last major or minor gas-rich interaction, then we 
would expect a rather tight correlation between age and asymmetry. We suggest that the source of extra scatter 
is natural diversity in progenitor star formation history, gas content, and merger mass ratio, but quantitative 
confirmation of this conjecture will require sophisticated modeling. We conclude that the asymmetry-age 
correlation is in basic accord with the merger hypothesis, and indicates that an important fraction of the early- 
type galaxy population is affected by major or minor mergers at cosmologically-recent times. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — galaxies: interactions — 
galaxies: peculiar galaxies: evolution — galaxies: general 
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1. INTRODUCTION 

It has become clear in the last decade that the distribution 
of galaxies is bimodal in color, and that this bimodality cor- 
relates with the structu re (or morphological type) of galax- 
ies (IStrateva et al.ll2001l:lKauffmann et al.ll2003l:lBlanton et al.1 
|2003[) . Early-type galaxies are dominated by a spheroid 
of stars supported largely by random motions, and tend 
to have red colors and very low present-day star forma- 
tion rates. They populate a tight color-magnitude rela- 
tion (CMR ; Sanda ge & Visvanathanlll978llBower et al.lll992L 
iRuhlandet al.1 120091) : the tightness of this CM R is an indi- 
cator of their rather ancie nt stellar populations dTrager et alJ 
2000; Gallazzi et al. 2006). In contrast, late-type galaxies are 
dominated by a rotationally- supported disk of stars, and tend 
to have significant ongoing star formation at the present day, 
and populate a reasonably broad (at optical wavelengths) blue 
cloud whose color s reflect a combination of recent star for- 
mation and dust (ISchiminovich et al1l2007|) . The stellar mass 
functions for the two classes of galaxy are different; early- 
type galaxies can extend to mu ch higher stellar mass than late- 
type galaxies (Bell et al. 200 3|). This b im odality persists at all 
redshifts up to z - 1 dBell et al.ll2004allbl : iFaber et al.ll2007h . 

One of the most remarkable features of the bimodality in 
the galaxy population is the observed correspondence be- 
tween galaxy structure and stellar populations. Models of 
galaxy formation in a cosmological context predict that the 
structure of a galaxy re flects the dynamical assembly (merger 
history) of a ga laxy dkauffman n. White & Guide ronl 11993b 
ICole et afll2000|) . whereas the cooling of gas and subsequent 



star formation is a property of the warm/hot reservoir of gas 
in a h alo, which is rel a tively unaffected by galaxy interac- 
tions (ICole et al] 120001: iCattaneo etaLl [2006). Galaxy for- 
mation models with 'standard' ingredients — the formation 
and evolution of dark matter halos, gas cooling, star forma- 
tion and stellar feedback — predict that almost a ll galaxies 
shoul d be forming stars at an appreciab le rate (ICole et al] 
120001: iCroton et all 120061: ICattaneo etaLl 120061: iKeres et al.1 
120091) . This is in direct disagreement with observations of 
a red sequence, and is usually interpreted (for good reasons) 
as indicating not just a minor failing of (admittedly uncertain) 
prescriptions for these physical processes; rather, this defi- 
ciency is interpreted as the signature of a completely separate 
physical mechanism. 

In this regard, the empirical correlation between a lack 
of star formation and a prominent stellar spheroid (bulge) 
is an imp ortant clue to the m echanism suppressing star 
formation ( Kauffmann et al. 2006). A bulge appears to be a 
necessary (but no t sufficient) requirement for the shut-down 
of star formation (Bell 2008). Taken together with the bulge 
mass-black hole m ass correlation (Magorrian et al. 1998; 
Hari ng & Rixll2004h . it is not unreasonable to postulate that 
feedback from accretion onto a supermassive black hole 
(active galactic nucleus [AGN] feedback) is an important 
mechanism by w h ich ga la xies quench their star formation 
dDi Matteo et al.1 [20051; iSpringeL Pi Matteo & Hernauistl 
1 2005 1 : ICroton et al.1 [20061: iBower et a l. 2006; Cattaneo etall 
pOOol: iHopkins et al.1 I2008allbl : iSomerville et al.n 2008: 
iJohansson. Naab & Bur kert 2009|); although oth er physical 
mechanisms may also play a role (see, e.g., Naa b et al.l 
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2007, Dekel & Birnboim 2008, Khochfar & Ostriker 2008, 
Johansson, Naab & Ostriker 2009, Dekel & Birnboim 2006, 
Birnboim et al.l 120071 or iGuo & Qhl 120081 on gravitational 



heating, the influence of the development of virial shocks, 
and the heating of large halos with cosmic ray energy). 

Galaxy merging play s a decisive role in this pic ture. 
iToomre & Toomrd (119721) and lBarnes & Hernquistl (119921) ar- 
gued that a natural formation route for early-type galax- 
ies (in particular elliptical galaxies) was through the ma- 
jor merger of two pre- e xisting galaxies (altho ugh see, 
e.g., lBurkertetal.1 120081 iNaab & Ostrikerl 120091 for dis- 
cussions of limitations of this obviously over-simplified 
picture). In the intervening time, an impressive vari- 
ety of evidence has been amassed showing that the prop- 
erties of early-type galaxies, at least at the broad level, 
are consistent with thos e expected for merger remnants : 
e.g., overall kinematics ([N aab, Khochfar, & Burkert 2006; 
Uesseit et al .1120071: Hoffman et al. 20Q|), the widespread ex- 
istenc e of kinematically-decoupled cores (Emsellem et al.l 



2007), the distribution of isophota l shape s (INaab & Burkertl 
20031: INaab. Khochfar. & Burkertl 120061). surface bright- 



ness p rofile shapes ([Naab & TrujilloT T2006l: IHopkins et al.l 
2008q), and low-level tidal debris around many elliptical 
galaxies dMal in & Carter 1983; Sc hweizer & Seitzed fl99l 
Ivan Dokkuml 12005: Tal et al.ll2QQ9h . Observations of merger 
remnants support this picture also, having surface brightness 
profiles, velocity dispersions and sizes in the near-infrared 
that are similar to intermediate-luminosity early- type galax- 
ies (iDasyra et al.l l2006: Roth eberg & Josephll2006|) . 

It is argued that shocks and non-circular motions in a galaxy 
merger cause the gas to lose angular momentum, leading to 
the growth of the black hole mass and conditions suitable for 
efficient c oupling of energy from an AGN to the diffuse gas in 
a galaxy (ISanders & Mirabellll996l: IHopkins et al.ll2008allbh . 
Circumstantial support for this picture of merger-driven AGN 
activity and feedback is seen in the tendency for the AGN 
hosts to h aye intermediate colors (betwee n blue cloud and red 
sequence; ISchawinski et al.ll2007l |2010|) . rapidly-outflowing 
metal-enr iched gas at ~ lOOkpc fr om an active quasar at 
z ~ 2.4 dProchaska & Henna wi 2009), and rapidly-outflowing 
gas (> 1000 km s -1 ) in post-starburst galaxies — such veloc- 
ities are not expected to result from starburst-driven winds 
and were interpreted a s being relics of a quasar-driven wind 
dTremonti et al.l 120071) . In such a picture, the suppres- 
sion (or quenching) of star formation is reasonably rapid 
and follows the merg er event ([Kauffmann & Haehnelt 2000; 
iDi Matteo etaI1l2005l) . 

Another aspect that has come into focus in the last years 
is that the formation of red sequence (and early-type) galax- 
ies is an ongoing process; in particular, more than half of the 
early-type galaxy population has come into place since z 



1 (se e e.g.. iBell et al.ll2004aL iBell et al.ll2004bl iFaber et al.l 
I2007L iBrown et al.l 12007b see |yan Dokk um et alJ 120 101 and 
iRobaina et al.ll20ldl for analyses focused on the most massive 
galaxies). Thus, one expects that the mergers (that create the 
early-type galaxies) and the physical processes that lead to 
the truncation of star formation (that make those early-types 
red) continue to happen at relatively recent times. In this 
case, these physical processes, and their late-ti me signatures, 
should be observable at the present day (see IHopkins et al I 
l2010al for a detailed and thoughtful discussion of merging and 
early-type galaxy evolution in a cosmological context). 
In this context, the work of Schweizer & Seitzer (1992) is 



of much importance. Motivated by the merger hypothesis for 
the origin of early-type galaxies, they sought a correlation be- 
tween fine structur e (from a tidal orig in) and color offset from 
the CMR (in Schweize r et al.l [T990 they did a similar work 
with a smaller sample with absorption line indices). If gas- 
rich mergers are the primary way in which early-type galax- 
ies are made, then one e xpects a correlation betwee n color 
offset and fine structure. ISchweize~ & Seitzerl (Il992|) found 
such a correla tion with ~ 4a c onfidence, albeit with consid- 
erable scatter; iTal et al.l (|2009|) found a similar result with a 
compl ete sample of nearby ellipticals. ISchweizer & Seitzerl 
(1992) used stellar population models to explore this result in 
more detail, finding that the early-type galaxy colors could be 
reproduced by a variety of post-merger ages from 2 to 6 Gyr, 
depending on the model parameters used. 

Yet, a number of issues are left open by these works. 
Number statistics was a clear issue in both cases (e.g., 
ISchweizer & Seitzerl 1 19921 had a total of 69 galaxies in their 
sample; ITal et al.l 120091 had 55 galaxies in their sample), 
and testing with a larger sample is cl early desirable. The 
definit ion of fine structure adopted by ISchweizer & Seitzed 
(1992) is motivated by the type of ti dal features regularly 
seen around nearby e a rly-typ e galaxies (jMalin & Carter 19831: 
Sch weizer & Seitzed [19881: Ivan Dokkuml 120051 : iTaletalJ 
2009), but has two limitations. Firstly, it is imp ossible to 
meas ure their fine structure parameter automatically. ITal et al.l 
( 2009) partially remedied this by defining a reproducible met- 
ric that is readily applicable to deep data, although they stud- 
ied a sample of elliptical galaxies only. Secondly, high- 
contrast tidal features are often the result of minor mergers 
or accretions, not major me rgers. Finally, at the time at which 
ISchweizer & Seitzerl (119921) was written, it was impossible to 
disentangle age and metallicity influences on spectral indices 
or colors; it is important to verify that age is in fact the driving 
parameter of any correlations seen. 

Cognizant of these issues, we have initiated an ef- 
fort to confirm and expa nd on the landmark results of 
ISchweizer & Seitzerl (119921) . We have defined a volume and 
^-band luminosity limited sample (® of 0.01 < z < 0.03 
early-type galaxies using the T wo Micron All- Sky Survey 
(2MASS; Skrutskie et al. 2006), the Sloan Digital Sky Sur- 
vey (SPSS: lAbazajian et al.ll2005|) . supplemented with red- 
shifts from the NASA/IPAC Extragalactic Database (NED). 
We have used asymmetry as an automated and repeatable met- 
ric for tidally-induced structure, and have analyzed it in con- 
junction with SDSS colors, and, crucially, ages and metallici- 
ties d erived from absorption-line spectroscopy ( Ga llazzi et al.l 
120051) : these parameters are described in ^3j Because of our 
use of the SDSS as our imaging data we analyze the brighter 
inner parts of early-type galaxies (less sensitive diagnostic of 
faint tidal features) but with the advantage of a larger sample 
(over 600 galaxies). We present our results in ^4j Finally, 
bearing in mind the insight gleaned from redshift surveys of 
red sequence and early-type galaxy evolution, we were in a 
position to model the spectr al properties of early-type galax- 
ies more completely than Sch weizer & Seitzerl (119921) : this 
discussion is presented in §0 The casual reader is invited 
to skip to §E1 directly, focusing in particular on Figs. [6] and 
ITTI In what follows, we use AB magnitudes and assume 
H = 70kms- 1 Mpc" 1 , ft m , = 0.3 and ^ A , = 0.7. 

2. DATA AND FITTING 
2.1. Sample selection 
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Fig. 1. — The estimated completeness of a 2MASS K-bmd selected sam- 
ple, as assessed by comparison with synthesized K-bmd apparent magnitudes 
of a deeper optically-selected catalog from the SDSS. The 2MASS K-selected 
sample is less complete for blue objects. 
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Fig. 2. — K-band magnitude dependence of redshift completeness of the 
sample. The three histograms show the redshift completeness of the 2MASS- 
selected sample using data from different sources: The dotted line shows 
NED completeness only, the dashed line SDSS completeness, and the solid 
line shows NED and SDSS combined completeness. 



Our aim is to test the relationship between fine structure 
and color deviation from the CMR for early-type galaxies. In 
order to test this relationship, a representative (but not nec- 
essarily complete) sample is required, and at distances close 
enough to allow us to discern both the global and tidal struc- 
ture of early-type galaxies. Accordingly, we have chosen 
galaxies within a relatively thin, and nearby, distance shell 
0.01 < z < 0.03 for study. The bulk of our sample can be 
drawn from the SDSS, but the relatively bright galaxies in this 
sample lack SDSS spectra (for a variety of reasons, but par- 
tially because their fiber magnitudes were too bright to target 
as part of the SDSS main survey). 

Accordingly, we have taken a hybrid approach for this 
work. The primary object selection was made from the 
2MASS All-Sky Extended Source Catalog. We selected 
galaxies brighter than 13.5 mag in K. We retrieved redshifts 
for these objects either from NEE0 or SDSS, matching the 
objects by position; NED was used to fill in objects lacking 
spectroscopy from the SDSS. In order to achieve a good red- 
shift completeness we only included galaxies from regions of 
the sky that are fully covered by spectroscopic sample of the 
SDSSDR3. 

As described in detail by iMcIntosh et al.l (|2006|) . from a 
comparison to the SDSS main galaxy sample, it turns out that 
the K-band incompleteness of 2MASS is modest and the only 
major source of incompleteness are blue low surface bright- 
ness disk galaxies. In order to check the completeness of our 
pa rticular 2MASS selection we follow a simil ar philosophy 
to iBell et alJ (12003b and IMcIntosh et ail (120061) . by explor- 
ing the 2MASS properties of galaxies selected in the optical 
regime. We selected a sample of galaxies with known red- 
shifts and SDSS photometr y from NED . With the same pro- 
cedure as described in lBell et all (|2003|) , the SDSS u,g : r,i : z 
magnitudes and the redshifts were fit with stellar population 
template s pectra derived using the PEGAS E stellar population 
code fsee lFioc & Rocca-Vblm erange 19971 for an earlier ver- 
sion of the model). These template fits were used to estimate 
^-corrections for the sample, and permitted estimation of the 
K-band magnitude of a galaxy from the optical fluxes alone. 
In this way we could explore whether 2MASS detected all the 
galaxies which one would a priori expect to be detected. The 
integrated incompleteness is approximately 5%. Fig. [T] shows 

1 The NASA/IPAC Extragalactic Database, 

|http : //nedwww . ipac . caltech . edu/| 



the color dependence of the completeness at 0.01 < z < 0.03, 
showing that the sample is less complete at bluer colors. 

Figure [2] shows the completeness of the redshift data com- 
ing from the two sources. SDSS is most incomplete at the 
bright magnitudes, whereas in NED the completeness de- 
creases with K magnitude. This trend dominates the com- 
bined catalog, but the decrease of the completeness is much 
more moderate. With our selection criteria, at least one of 
the redshift data (SDSS or NED) is available for 91% of the 
K < 13.5 galaxies, the incompleteness mainly coming from 
faint, K > 12 galaxies. 

For the structural studies we select a nearby subsample 
which is volume limited between 0.01<z<0.03. The cor- 
responding and K-band absolute magnitude limit is M K < 
-22.31 (whic h approximately corresponds to L > L*/5, 
IBell et al.1l2003[) . This sample contains 3006 galaxies. Its esti- 
mated completeness is approximately 90%; this completeness 
is dominated by the spectroscopic incompleteness (largely 
from SDSS, and this is primarily a geometric effect owing 
to the inability to place fibers on all objects) coupled with a 
modest trend to miss out low surface brightness, blue objects; 
for red galaxies, the focus of this work, the sample is essen- 
tially limited only by SDSS fiber placement. 

2.2. Fitting 

We carried out image fitting with GALFIT to select early- 
type candidate systems. For fitting the galaxy profiles we con- 
structed an automa ted pipeline involving programs SExtrac- 
tor version 2.3.2 jB ertin & Arnoutslfl996h and GALFIT ver- 
sion 2.0.3. (iPeng et al.ll2002h . SExtractor creates a catalog of 
sources from an astronomical image. The output parameters 
can be used to create an initial setup for light profile fitting 
with GALFIT. 

We pipeline was run on SDSS r-band images. For each im- 
age, SExtractor was used to look up all sources and to identify 
the target galaxy, using SDSS positions as input. Besides our 
primary galaxy, we fit the nearby objects that could influence 
the light profile fit of the main galaxy. The selection of these 
secondary targets was made by distance and intensity: the ob- 
ject was ignored if the distance between the two objects ex- 
ceeded 1.5 times the sum of their sizes, or the intensity of the 
secondary object was lower than 1% of the main object. In 
some cases SExtractor identifies parts of the galaxies, such as 
large Hii regions as separate objects. In order to correct these 
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cases, check images were created with labeled fit targets. The 
images were visually inspected and unnecessary sources were 
untargeted. These results were used to create the GALFIT ini- 
tial setup. 

The main target galaxies were fitted with a ISersid (119681) 
profile 

E,(r) = E^" c ^ 1/n - 1 (1) 

where r e is the effective radius of the galaxy, S e is the surface 
brightness at this radius, n is the Sersic index (in the special 
case when n = 1 we get the exponential profile, n = 4 corre- 
sponds to de the Vaucouleurs profile) and c = l.9992n~ 03271 
so that the half of the flux is within r e . The nearby objects 
(secondary targets) were fitted either by Sersic or, if more 
star-like, by the point spread function (PSF). The PSF was 
created from stars of each field image separately, based on 
Moffat profile fits in the form: 

S ^ r) = r7T7~7 ^2^m ' ( 2) 

where is the dispersion radius and m is the power-law in- 
dex. 

Initial values for sky level, flux, size and axis ratio were set 
using SExtractor output parameters. The Sersic index and axis 
ratio for galaxies were set to initial values common for all ob- 
jects, n = 1.5, b/a = 0.9. The fit parameters were constrained 
using empirical (observational) and SExtractor output infor- 
mation: n was required to lie in the interval 0.5 < n < 7, axis 
ratio 0.3 < b/a < 1, position within 2 pixels of the SExtractor 
value, and half-light radius between and 1.5rsExtractor- The 
pipeline first fits the objects separately, keeping the other ob- 
jects fixed, and then finally fits the objects all together in the 
end, using their separate fit results as the initial guess. The 
sky is fixed at all times. 

One way to test the fitting pipeline is to fit simulated objects 
of known parameters. We used a suite of simulated galaxies 
composed from a Sersic bulge + exponential disk, with a wide 
range of magnitudes, bulge to total ratios and light profile pa- 
rameters. We selected a test sample similar to our sample of 
interest: bulge dominated systems, having bulge to total ratios 
B/T > 0.6 and magnitudes brighter than 16.5 (the K < 13.5 
cut corresponds approximately to r < 16.4). The simulated 
galaxies were fit by the pipeline described above. The fits give 
us an important sanity check of the fitting procedure. Our re- 
sults show that Sersic fit works relatively well, the simulated 
magnitudes are reproduced with an error ~ 0.2 mag for the 
test sample. 

When fitting the real images, the Sersic fit is successful and 
gives results similar to the SDSS photometric parameters for 
the majority of our galaxies. It also corrects the known error 
of the SDSS imaging pipeline that it tends to split up large 
objects into multiple sources. This affects 4.6% of the objects 
in our sample, mostly objects with n < 1 . The Sersic fit fails in 
1 .5% of the cases because the fit did not converge or the image 
file was corrupt. 1.9% of the fit results are unreliable because 
the fit did not move out from the initial setup or reached the 
fitting boundary. Figure [3] shows the fitted sizes (circularized 
r° = r e yjbjd) compared with the SDSS radii. The fit gives 
results consistent with the SDSS sizes within 15%, when > 3a 
outliers are discarded. 

We choose n > 2.5 galaxies for further study, which ap- 
pears to be a reasonable way to automatically rej ect a sig- 
nificant fraction of the disk-dominated galaxies (Bel l et all 
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Fig. 3. — The relative difference of the fitted circularized Sersic r e vs. 
SDSS Rsq as a function of r magnitude. Small black symbols show n < 2.5 
galaxies (a proxy for disk-dominated galaxies; n = 1.5 excluded), while small 
gray symbols show systems with n > 2.5 (a proxy for bulge-dominated galax- 
ies). Larger open circles show systems where the fit hit the upper limit in n, 
and crosses show where the fit hit the lower limit in r e . 

2004b fl We further exclude all objects with failed fits, as de- 
scribed above. The criterion n > 2.5 is deliberately generous, 
and includes a larg e number of galaxies with prominent disks 
dBell et al.ll2004bh . Accordingly, we classified the n > 2.5 
subsample by eye, classifying 553 objects as elliptical (pure 
bulge; E), 498 as lenticular (smooth, non- star-forming disk; 
SO), 479 as later types (Sa-Irr) and 27 as mergers. 

We choose to visually exclude all galaxies with stars or 
companion/projected galaxies which made a major contribu- 
tion to the asymmetry in the area in which asymmetry is mea- 
sured. The spirit of this cut is to excise all systems where we 
were concerned that the asymmetry measurement would be 
corrupted by a non-interacting (or non- strongly-interacting) 
object (be it a projection or an interaction). While some 
form of masking might have been desirable, we adopted this 
conservative approach to account for objects missed by SEx- 
tractor and because the extent of the masked area is never 
straightforward to choose. We lose a little over 1/3 of the 
sample using this cut (it is our most important cut), and the 
loss is primarily geometric (and therefore introduces no bias) 
for the galaxies lost from bright stars and projections, and 
those galaxies thrown out because of a pre-merger satellite. 
We have confirmed that the color-magnitude diagram of these 
excised galaxies is very similar to the final sample of galaxies 
without projections (there are a few more outliers in the red- 
wards and bluewards directions in the excised sample com- 
pared to the final sample, as is expected from a sample of 
galaxies where there is a projected star or galaxy contaminat- 
ing the photometry). There is some natural fuzziness in this 
process, and some advanced mergers where the second nu- 
clei (and large amounts of tidal disturbance) were within the 
asymmetry measurement area were included in the sample. 
After this final selection our sample contains 288 E and 319 
SO type galaxies, together 607 objects. 

3. ASYMMETRY AND STELLAR POPULATION MEASUREMENTS 

For the investigation of the correlation between the post- 
merger structure and stellar population we need to construct 
two broad classes of metric: descriptions of the (hopefully 

2 If instead a selection by r-band concentration as calculated by the SDSS 
(defined as the ratio of the aperture containing 90% of the flux in t he Petrosian 
radius to the half-light radius) being larger than 2.6 ( Strateva et al. 200l|) was 
used, the sample would remain unchanged at the 90% level. 
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Fig. 4. — Two examples of asymmetry measurement. Left panels: the orig- 
inal image, with the elliptical annuli over-plotted; right panels: asymmetry 
images. 



tidally-induced) non-equilibrium structure of an early-type 
galaxy ( 33. ]} , and parameters character izing the stellar popu- 
lation, like color, age, and metallicity ( 33.21) . 

3.1. Structure 

We choose to use asymmetry as a repeatable diagnostic of 
non-equilibrium structure. This measure is insensitive to sym- 
metric bar s and rings (u nlike the residual structure measures 
used by iTal et al.l 120091) . which are important structures in 
some SOs but are not (at least unique) indicators of a recent 
perturbation. Our metr ic is not as tailored a s the fi ne struc- 
ture parameter used by [Schweizer & Seitzer] (11992b , but the 
asymmetry measurement is repeatable and can be automati- 
cally applied to survey datasets. 

We ado pt an asymmetry measure similar to m ethods dis- 
cussed by IConselice. Bershady. & Jangrenl (120001) . We mea- 
sure asymmetry within an elliptical annulus of 0.8-3 r e . In 
this way we minimize the impact of the bright galaxy core 
on our asymmetry measure which would otherwise be very 
sensitive to centering uncertainties and the influence of dust 
disks and features in the inner part of early-type galaxiefl A 
different choice of the aperture does not influence the results 
qualitatively, unless the inner 0.5r e is included in the asymme- 
try measure. The asymmetry calculating algorithm is carried 
out as follows in the steps 1-7, for an illustration see Figure[4] 

1. Rotate the image around the GALFIT center by 180°. 

2. Subtract the rotated image from the original one. 

3. Integrate flux in quadrature within the elliptical annu- 
lus between 0.8-3 r e . This measure will have an ad- 
ditive contribution from noise which can be separated 

3 Note that initially we chose to adjust the center to minimize asymme- 
try within 0.5r e . Galaxies with significant asymmetry often showed that the 
center of the brightest isophotes was offset from that of the fainter isophotes, 
yielding a change in the asymmetry with radius. We found that the GALFIT 
center was very close to a choice that minimized A, with good S/N. Accord- 
ingly, in what follows, we choose to use the GALFIT center, and omit the 
central regions from the asymmetry measure. 



and corrected for. Let us separate the flux / into signal 
and noise term, / = s + n and analogously for the ro- 
tated image f = s' +n' . Then the asymmetry flux /-/' 
integrated in quadrature over the image will read 

5>-/') 2 =I>+«-*'-«') 2 

=E^') 2+ E" 2+ E« ,2 > ( 3 ) 

as signal and noise, as well as the both noise terms are 
uncorrected. This will result in a non-zero asymmetry 
even in the case of s = s', i.e. perfectly symmetric un- 
derlying signal. We can correct for this by estimating 
the noise contribution and subtracting it from the inte- 
gral. We estimate n-n' using the GALFIT model image 
with noise similar to the observed data. 

4. Apply Poisson noise + read noise to the (symmetric) 
GALFIT model image. 

5. Measure asymmetry in the same way as for the actual 
image (steps 1-3, rotate, subtract, integrate in quadra- 
ture). 

6. Subtract model asymmetry from the object asymmetry. 

7. In the end we make the asymmetry measure indepen- 
dent of the total flux: normalize by the object (sky sub- 
tracted image) flux squared and integrated within the 
same area. 

To summarize, our asymmetry measure is defined as 

a = (Ec/-/') 2 -^™-™') 2 ) /E</ 2 >- w 

where / denotes the flux of the galaxy, m is the fitted Ser- 
sic model, the primes denote the corresponding quantities ro- 
tated by 180 degrees, and the summation goes over pixels 
within the elliptical annulus 0.8-3 r e . Because the asymmetry 
is determined from the square of the asymmetry image, val- 
ues of ~ 0.01 are already very highly significant, suggesting 
asymmetries involving ~ 5% of the flux (this definition counts 
asymmetric flux from both the positive and negative parts of 
the asymmetry image; other definitions divide the asymmetry 
by two to account for this effect). 

The uncertainties in asymmetry were calculated in two 
complementary ways. Random uncertainties were estimated 
by subjecting symmetric model galaxies with realistic model 
and sky noise properties to the same procedure above; these 
random uncertainties are 5 A ~ 0.0003, and are shown by the 
horizontal component of the cross-shaped error bars in Figs. 
[6l[TT]and[T2l The second estimate of asymmetry uncertainty 
is systematic in nature, reflecting our choice to neglect the 
inner 0.8 r e when estimating the asymmetry (in order to min- 
imize the effect of small centering uncertainties, and the in- 
fluence of small dust lanes or disks in the inner parts of these 
early-type galaxies). We changed the inner cutoff to 0.6-1 r e 
and the outer cutoff to 2-4 r e ; the result was that the asymme- 
tries changed significantly, because of the different area being 
probed, but to first order all galaxies moved to higher or lower 
values of asymmetry together. Experimentation revealed that 
the best way to visualize this source of uncertainty is as an 
uncertainty in asymmetry rank. The rank changed typically 
by ~ 5%; we assign an asymmetry rank error of 5%, and 
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for the purposes of visualization we translate this rank uncer- 
tainty into an asymmetry uncertainty using the asymmetries 
of galaxies that have ranks 5% different from the galaxy of 
interest. We show error bars calculated in this way in Figs.[6l 
[TT] and El as the solitary horizontal error bars. The A uncer- 
tainties from this source are a strong function of A, increasing 
from very small values for low values of A to large values for 
the tail of galaxies extending to large A. 

3.2. Stellar population diagnostics 

We choose for consideration in this work two types of 
stellar population diagnostics: galaxy colors, as an in- 
dicator of broad changes in the stellar population; and 
luminosity-weighted stellar ages and metallicities — avail- 
able for roughly 2/3 of the sample — as model-dependent but 
more information-rich stellar population diagnostics. 

We use SDSS model g-r color corrected for foreground 
extinction and A B offset, then k-corrected using method of 
iBell et al.1 (|2003[) . Note that these colors are very similar 
to colors we measured ourselves within < r/r e < 3, and 
the results change only imperceptibly when such colors are 
used. We adopt the model colors because these are read- 
ily and publicly available. We assign uncertainties of 0.03 
mag in g-r color to account for random photometric errors 
(sub-dominant) and calibration uncertainties from galaxy-to- 
galaxy in the zero-points of the SDSS photometry (the domi- 
nant contribution). The color combination g - r was used be- 
cause it is the bluest well-measured SDSS color combination. 
While, in principle, u — r color is more sensitive to stellar pop- 
ulation (as it straddles the 4000 A break), w-band has poor S/N, 
and suffers from increased systematics compared to g-band. 

In order to aid in disentangling the effects of the age/star 
formation history of a galaxy from metallicity in driving its 
color, we u se estimates of lumino sit y-weighted age and m etal- 
licity from lGallazzi et al.l (120051) . iGallazzi et all (120051) de- 
rived ages and metallicities from all SDSS spectra with me- 
dian S/N> 20 per pixel. This stringent S/N cut selects only 
~ 1 /4 of all galaxies, primarily those with high surface bright- 
ness cores. Our sample consists of E/S0 galaxies only, and 
therefore is much more complete than the whole sample: 97% 
of the Es and SOs in the sample with SDSS spectra i n DR 4 
have age and metallicity estimates from Gallaz ziet al.l (|2005). 

Ages and metallicities are derived by comparing observed 
spectral absorption features (three primarily age-sensitive in- 
dices, 4000 A break strength D4000, H/3, H 7a + H5 a , and two 
primarily metallici ty- sensitive indices [Mg2Fe] and [MgFe]') 
with a library of Bruzual & Chariot (120031) models with a 
wide range of star formation histories and metallicities. The 
models have solar abundance ratio patterns (i.e., they do not 
allow for a-element enriched abundance patterns), but the im- 
pact of this shortcoming on the derived ages and metallici- 
ties is not that severe owing to the choice of combined Mag- 
nesiu m and Iron indices as the metallicity- sensitive indices 
(see IGallazzi et al.|[2005l for details and more discussion). In 
this paper, the main emphasis is on the age estimates (re- 
call that the hypothesis being tested is that richly- structured 
early-type galaxies have younger ages), which depend pri- 
marily on the B aimer lines and the strength of the 4000 A- 
break. These ages are luminosity-weighted, i.e., the ages are 
the mean luminosity-weighted ages of the models that pro- 
vided a reasonable fit to the five observed spectral indices, and 
are correspondingly younger than the mass-weighted ages of 
those model populations. 



-23 -22.5 -22 -21.5 -21 -20.5 -20 -19.5 -19 

M r,model 

Fig. 5. — The CMR using g and r-band SDSS model magnitudes for galax- 
ies in our sample with measured asymmetries and meaningful SDSS model 
magnitudes. The line shows the best fit to these data. 



Section 2.4 of Gall azzi et al.l (I2005I) discusses in detail the 
sources of random and systematic error in these age and 
metallicity estimates. Briefly, random errors (from the spec- 
tra themselves) do not cause any systematic shifts in the me- 
dian ages or metallicities from their correct values (verified 
using simulations), and for our subsample lead to a typical 
uncertainty in log age of 0.1 dex and metallicity of 0.09 dex. 
The main systematic errors were identified to be alpha ele- 
ment overabundance uncertainties and the degree of bursti- 
ness in the star formation history of the galaxy population. 
The systematic errors in metallicity and log age from these 
sources was argued to be less than 0.05 dex for relatively large 
variations in alpha element overabundance and the fraction of 
galaxies undergo ing l arge bursts (we will discuss further these 
uncertainties in 34.11) . It is worth noting that the luminosity- 
weighted ages are reasonably model-dependent, and depend 
on when star formation is assumed to start and to the tim- 
ing and amplitudes of bursts of star formation. For example, 
a model with given line strengths could be fit by a popula- 
tion that started forming stars 12 Gyr ago, or only 5 Gyr ago; 
this will naturally lead to dramatic differences in luminosity- 
weighted age. Accordingly, in this work, we use luminosity- 
weighted ages primarily as a tool to understand relative trends 
in age, and to attempt to disentangle the effects of age and 
metallicity on broad-band colors. 

4. RESULTS 

Figure [5] shows the CMR for the 605 E/S0 galaxies in our 
sample with measured asymmetries and SDSS Model magni- 
tudes. We fitted a line to the main locus of the CMR, exclud- 
ing blue outliers at g — r < 0.6 from the fit. The equation of 
the fitted line is 



g-r = -0.023 M r + 0.28. 



(5) 



We subtract this linear relation from the color to get the color 
offset from the red sequence: A(g-r). The scatter of A(g-r) 
is 0.041 with 3a-clipping applied; much of this scatter is in- 
trinsic (not measurement error, the typical measurement er- 
ror in model g-r for our sample is ~ 0.03 mag and is dom- 
inated by flat fielding and calibration uncertainties; see also 
iSchweizer & Seitzedll992l and iRuhland et al. 2009 who also 
measure similar intrinsic scatter in the CMF0). The mean 



4 for the purposes of comparing scatter in ^ 
V)~2.5Afe-r) 



-r with U-V, note that A(U - 
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Fig. 6. — Color-asymmetry diagram. Left panel: Es in our sample, right panel: SOs in our sample. The li nes show averages in bins, with error bars denoting 
the error of the average. The cross-shape error bars indicate the random error of asymmetry at A=0 (see Sec. l3.U and the color error. The systematic uncertainty 
due to the choice of the asymmetry measurement aperture is shown as solitary horizontal error bars at two different asymmetry values. A correlation between 
color and asymmetry is present in both galaxy types. 



color offset is slightly negative, (A(g-r)) = -0.005, because 
of the clipping of very blue outliers from the CMR fit. 

Figure [6] shows the color offset A(g-r) as a function of 
our asymmetry measure, with bin averages plotted with lines 
(5 bins equally spaced between A of and 0.01, a bin with 
0.01 < A < 0.02 and one with 0.02 < A < 0.04). The data 
show a correlation between asymmetry and color offset for 
both E and SO galaxies, with approximately the same slope 
and zero point; this figure can be compared reasonab ly di- 
rectly with Figs. 2 and 3 of Sch weizer & Seitzerl {1992), and 
confirm their detection of a significant correlation between 
substructure in early-type galaxies and the offset from the 
CMR. 

One obvious feature of this correlation is that it is very scat- 
tered; the scatter in the CMR exceeds the magnitude of the 
correlation between color offset and asymmetry. This scatter 
is pri marily intrinsic (i.e., not measurement error; see for ex- 
ample Gallazzi et al. 2006 for a demonstration that the stellar 
populations of galaxies on the red side of the CMR differ sub- 
stantially f rom tho se on the blue side). We explore this issue 
in sections I47TI and [5X4l in some depth. 

The correlation is clearest for A < 0.01 — the run of both 
the mean and median color at A < 0.01 is very similar. At 
larger asymmetry values, in particular for elliptical galaxies, 
there are quite a few galaxies with significant asymmetry but 
relatively red colors — the median color for A > 0.02 el- 
lipticals is reasonably red. We discuss this later in 35.21 in 
the context of dissipationless merging between already gas- 
poor early-type galaxies. Nonetheless, the bulk of the early- 
type galaxy population (those with A < 0.01) shows a trend 
between color offset and asymmetry. To help visualize the 
trend, we ranked the E/S0 galaxies by A and show how rank 
A correlates with color in Fig. We have investigated the 
significance of the correlation between A(g-r) and A us- 
ing Pearson's correlation coefficient r and Spearman's rank 
correlation coefficient s (see Table Q] for the full breakdown 
by galaxy type and stellar population diagnostic^- The cor- 
relation between A(g-r) and A exists but has large scatter. 

5 We present one-sided probabilities as we test the significance of the ob- 
served (and expected) anticorrelation between stellar population diagnostics 
(color/age/metallicity) and asymmetry. 



Over the full range of color and asymmetry, the correlation 
between A(g-r) and A has less than a 0.0035% probability of 
being from chance alone as estimated using Pearson's corre- 
lation analysis (or a 0.015% probability of being from chance 
alone if one restricts attention to relatively low asymmetries 
and small color offsets). If instead one studies ranked color 
offsets and ranked asymmetries, a Spearman's rank correla- 
tion analysis argues that the relationship between asymmetry 
and color for the combined E/S0 sample has less than 0.5% 
probability of being given by chance alonfi In all cases, the 
values of correlation coefficient and significance change lit- 
tle if different contributions of inner/outer radii are used for 
defining the asymmetry values. 

In Figures [U and [TO] we show several examples of dif- 
ferent asymmetry values for both morphological types, to 
allow the reader to generate some intuition for the asym- 
metries shown by the sample. In Fig. [8] shows a sample 
most symmetric Es and SOs, with asymmetry in the range 
< A < 0.001. Fig.[9]shows galaxies with 0.005 < A < 0.01; 
galaxies in this range are those that drive the bulk of the 
correlation between color offset and asymmetry. Fig. [TO] 
shows highly asymmetric galaxies from the asymmetry bin 
0.02 < A < 0.04. It is clear that the asymmetries being dis- 
cussed here are not from isolated Hll regions, or low-level spi- 
ral arms; rather, these asymmetries are to a great extent driven 
by large-scale asymmetries, many of them matching qualita- 
tively the form expected for tidal features. Taken together 
with Fig. |6l these figures under line th e primary result of this 
work (and lSchweizer & Seitzedll992l:lTal et al.ll2Q09h : early- 
type galaxies with asymmetric signatures of tidal interactions 
tend to be offset from the locus of the CMR, in the sense that 
more asymmetric galaxies tend to be bluer than undisturbed 
galaxies (albeit with large scatter around this trend). 

4.1. What drives these correlations, age or metallicity? 

One of the unavoidab le limitations of the work of 
ISchweizer & SeitzeH (119921) was that it was impossible to ex- 
plore quantitatively whether it was age or metallicity that was 

6 This is a relatively modest significance, and if this were the only set of 
stellar population diagnostics, one may be concerned abo ut the existence of 
this correlation. We will illuminate this issue later in 34.11 
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TABLE 1 

Correlation Probabilities of stellar population-asymmetry correlations 



Sample N 


r P(t < fobs) 


s 


P(s < S bs) 


Color-asymmetry 


E 


288 


-0.203 0.00026 


-0.114 


0.027 


SO 


319 


-0.0851 0.065 


-0.111 


0.024 


all 


607 


-0.161 3.5xl0~ 5 


-0.111 


0.0031 




Color-asymmetry A < 0.01 and A(g-r)| 


<0.15 


E 


250 


-0.169 0.0037 


-0.107 


0.046 


SO 


275 


-0.139 0.010 


-0.121 


0.023 


all 


525 


-0.157 0.00015 


-0.114 


0.0045 


Age-asymmetry 


E 


190 


-0.37 7.4 xl0~ 8 


-0.363 


2.9 xl0~ 7 


SO 


199 


-0.316 2.8xl0~ 6 


-0.375 


6.4 xl0~ 8 


all 


389 


-0.351 5.4xl0- 13 


-0.376 


6.2 xlO- 14 


Metallicity-asymmetry 


E 


190 


-0.22 0.0012 


-0.0486 


0.25 


SO 


199 


-0.073 0.15 


-0.0178 


0.40 


all 


389 


-0.173 0.00031 


-0.0235 


0.32 



Note. — Probability of observing a given stellar population-asymmetry correlation in a completely uncorrelated dataset, as evaluated using Pearson's 
correlation coefficient r and Spearman's rank correlation coefficient (s). The probability of a given value of the correlation coefficient r being due to chance alone 
in an uncorrelated dataset is distributed as a Student's t distribution with N- 2 degrees of freedom, where t = rj \J (1 -r 2 )/(N -2). 




-50 50 100 150 200 250 300 
A rank 



-50 50 100 150 200 250 300 350 
A rank 



Fig. 7. — Color offset A(g-r) shown as a function of the rank of asymmetry. Left pan el: E s, right panel: SOs. Symbols are same as in Fig. [6] the cross shaped 
error bars show the typical systematic uncertainty of asymmetry (see explanation in Sec. I3.lt and the color error. 



driving the color-fine structure relation. Here, we make use 
of estimates of luminosity-weighted stellar age and metallic- 
ity, available f or almost all of the y 2/3 of our sample with 
SDSS spectra (Gallazzi et al.ll2005|) . to address this question 
quantitatively. 

Figs. Qj] and [l2l show the run of luminosity-weighted age 
and metallicity as a function of asymmetry separately for 
E and SO galaxies (the trends in median age/metallicity are 
very similar). There is a very strong correlation between 
age and asymmetry for both E and SO galaxies (combined 
E/S0 sample has an essentially zero chance of being due to 
chance alone, from a Pearson's or Spearman's correlation co- 
efficient), and a substantially weaker or no correlation for 
metallicity (with a 0.031% chance of being from chance alone 
according to a Pearson's correlation coefficient, and a 1/3 
chance of being from chance alone as measured by a Spear- 
man's coefficient). 

It is worth discussing briefly the sources of uncertainty in 



this result. The random uncertainties (~ 0.1 dex in both log 
age and metallicity) are of little importance in this particular 
analysis, as one is seeking relative trends within a large and 
homogeneous dataset. The composite datapoints (mean for 
subsamples of 20 or more galaxies, typically) have small ran- 
dom uncertainties (< 0.02 dex), and differences of ~ 0. 15 dex 
in e.g., age between A ~ and A ~ 0.01 are detected with 
great significance. Systematic effects could be more impor- 
tant. For example, if there were systematic differences in 
alpha element overabundance as a function of asymmetry to 
the tune of ~ 0.3 dex in alpha element overabundance (a very 
large effect), a trend in age of up to 0.05 dex could be pro- 
duced (smaller than the observed trend in age with asymme- 
try) in conjunction with an offset in metallicity of simi lar am- 
plitude but opposite sign (§2.4 of Gallazz iet al.ll2005l) . Such 
a tendency is not seen; it is likely that alpha element over- 
abundance is not an important driver of the observed trend in 
age (and weak or no trend in metallicity). Furthermore, as 
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Fig. 8. — Galaxy images (first and third column) with their corresponding 
asymmetry images in elliptical annuli (second and fourth column, respec- 
tively) for the < A < 0.001 asymmetry bin (most symmetric objects). First 
and second row: Es, third and fourth row: SOs. Each box is approximately 
55 arcsec on a side. A 10 kpc horizontal line is shown for each object. 




FIG. 9. — Same as Fig.[8]but asymmetry range 0.005 < A < 0.01. These 
are typical asymmetry values for the galaxies that are driving the color- 
asymmetry correlation. 



we saw earlier, changes in the inner/outer radii used to de- 
fine asymmetry also change the results very little. One is left 
to conclude that genuine differences in star formation history, 
either through offsets in the age of the bulk of stars in the 
systems, or through a higher incidence of significant recent 
(last few Gyr) star formation activity between symmetric and 
asymmetric early-type galaxies is the main driver of the age- 
asymmetry (and therefore color-asymmetry) trend. 

We can therefore conclude that the color-structure corre- 
lation is, to the extent that we can tell, driven primarily by 
age effects; thi s observation verifie s the fu ndamental assump- 
tion made by Schweizer & Seitzen (11992 ). The correlation is 
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Fig. 10. — Same as Fig.[8]but asymmetry range 0.02 < A < 0.04 - most 
asymmetric E/S0 galaxies. 



actually substantially clearer in age-asymmetry space, as the 
contribution from random scatter in the metallicities to the 
color-asymmetry trend act to scatter out the color-asymmetry 
relation. 

5. DISCUSSION 

Motivated by the work of ISchweizer & Seitzert (Q992), we 
showed that the color offset from the early-type galaxy CMR 
correlates with (primarily tidally-induced) asymmetry, al- 
though the scatter in the correlation is large. We showed, for 
the first time, that this correlation is driven primarily by age 
effects. In this section, we will explore further the meaning 
and limitations of these results. We will present a quantita- 
tive model-based discussion of the results, focusing on the 
bearing that these results have on the hypothesis that one of 
the primary formation routes of early-type galaxies is through 
galaxy merging. Then, we will wrap up with a discussion of 
the limitations of this work, and scope for future improve- 
ments. 

5.1. A quantitative model of red sequence offsets 

Like ISchweizer & SeitzeH (Q~992), we show a relationship 
between color offset fro m the red sequence and as ymmetric 
structure of a galaxy. In Schwe izer & SeitzeH (Il992|) . it was 
argued that the 'fine structure' being measured was a strong 
indication of a previous gas-rich galaxy merger. Thus, in that 
paper, they presented a model in which early-type galaxies 
were the product of gas-rich major merging, and they stud- 
ied the time taken for major merger remnants to age as a way 
of translating their color offsets into timescales for early-type 
galaxy evolution. Their model had a number of free parame- 
ters: the two key parameters and e controlled the duration 
and prominence of star formation after the merger event. They 
found that these parameters were of great importance in set- 
ting the timescales for fading of the remnant. In models with 
a large value of the remnant reddened only slowly, reaching 
the red sequence only after many Gyr. In models with short- 
lived post-merger star formation, they found a rapid evolution 
onto the red sequence. 

One of the motivations of their work was to address the 
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Fig. 11. — The logarithm of the luminosity weighted age (in years) versus asymmetry. Left panel: elliptical galaxies, right panel: lenticular (SO) galaxies. 



Symbols, lines and error bars are as in 
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Fig. 12. — Metallicity log lo (Z/Z0) as a function of asymmetry. Left panel: Es, right panel: SOs. Symbols, lines and error bars are as in Fig. [6] 



'gap' between clear remnants from relatively recent merg- 
ers (< 1 Gyr ago, e.g., NG C 3921 or NGC 7252) an d 'an- 
cient' early-type galaxies. ISchweizer & Seitzerl (Il992|) quote 
King's question to Toomre at the 1977 Yale conference!]: 
"You showed us 10 merging pairs and then asked us to look 
for, or at least accept the existence of, 500 remnants from so 
long ago that they no longer bear the 'made by Toomre' label. 
I would be much more impressed if you showed us the 20 or 
30 such systems in the box immediately adjacent in your his- 
togram. What do these mer ged pairs look like in their n ext 
few galactic years?". While lSchweizer & Seitzerl (119921) ar- 
gued that the blue, structured early-type galaxies in their study 
were reasonable candidates to fill the 'King gap', the consid- 
erable sensitivity of the evolution of the merger remnant to the 
properties of the merger and star formation after the merger 
made it difficult to reach a unique conclusion. 

In this subsection, we revisit this issue with the benefit of 
more advanced stellar population models, and more critically, 
a more physically-motivated model for the color evolution of 
early-type galaxies. We will argue that the systems explored 

7 Neither of us were at the Yale conference. 



in this paper (and in ISchweizer & Sei tzer 1992) fill the 'King 
gap', but that owing to natural scatter in the merger history 
of early-type galaxies and in the stellar populations of these 
newly-formed early-type galaxies, it is impossible to connect 
systems one-to-one with a given merger time with data of the 
type that we use here; King's question was ill-posed, and will 
be very difficult to answer unambiguously. 

5.1.1. Model ingredients 

Clearly, there is a huge parameter space that can be ex- 
plored when modeling the evolving early-type galaxy popu- 
lation, and arriving at a single unique picture that is clearly 
superior to the others will be impossible (as clea rly articulated 
and demonstrated by ISchweizer & Seitzer 1993). 

We take an approach that is motivated by the last decade 
of research on the evolution of the early-type galaxy popu- 
lation. There are now a number of surveys that have empiri- 
cally tracked the evolution of the stellar mass density in early- 
type galaxies, with the general result that the growth in stellar 
mass density is driven pri marily by an increase in the num- 
ber of early-type galaxies dC hen et al. 120031: iBdl et al.ll2004al: 
iBrown et a l. 2007: iFaber eT al. 2007). Star formation is shut 
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off in blue star-forming galaxies , and the remnant fades and 
reddens onto the red sequence (Bell et al. 2007). The exact 
mechanisms that shut off star formation are still a matter of 
some debate (merging and AGN feedback, environment, or 
gas exhaustion are all possibilities; furthermore, the balance 
of such processes and timescales on which they operate is 
likely to depend on a variety of parameters such as galaxy 
mass, or merger mass ratio; e.g., Johanss on. Naab & Burkert 
2009); also open is the issue of how important any burst of 
star form ation might be bef ore the truncation of star forma- 
tion (see Rob aina et al.l 120091 for a discussion of the average 
effects of galaxy mergers on star formation rate). Notwith- 
standing this uncertainty about mechanisms for shutting off 
star formation, it is not unfair to assume that the timescale 
for truncation of star formation is short, <C 1 Gyfl Evi- 
dence indicating the short timescale for star formation trun- 
cation include the detection of ra pid, arguably AGN-d riven 
winds in post-starbursi ^galaxies (|Tremonti et al.l |2007|) . the 
pro nounced bimodalit> Bin the star formation rates of galax- 



ies (IStrateva et al.l 200 1 



Bell et al. 2004a; Schiminovich et al. 



|2007|) . the high incidence of concentrated, post-starburst AGN 
galaxies in the Val ley' between the star-fo rming and non- star- 
forming galaxies (|Schawinski et al.l2007|) . the modest popula- 
tion of blue morphologically early-type galaxies (if extended 
periods of star formation followed a merger, a large fraction 
of star-forming early-type galaxi es would be observed; such 
a pop ulation is not common; H auBleri I2007L iRuhland et al.l 
2009, this paper) and analys is of the spectra of early-type 
galaxies (Gall azzi et al.ll2006|) . In the context of this paper 
(which was motivated originally by the merger hypothesis), 
we will suppose that merging rearranges the stellar content of 
the remnant into a more spheroidal configuration and leads to 
the suppression of future star formation through AGN feed- 
back or some other mechanism. 

Thus, we model a constantly-growing early-t ype galaxy 
population (an exercise sim ilar in spirit to. e.g.. lHarker et all 
120061 IRuhland et all I2009L or iNaab & Qstriked 12009b asa 
population of galaxies with constant star formation rate, that 
then stop forming stars at times drawn from the following 
probability distribution: P(t) ex 1 -£(Gyr)/21, where t is the 
age of the galaxy in Gyr. For concreteness, we choose to 
assign to = 12 Gyr, i.e., the formation of early-type galaxies 
in this model started 12 Gyr ago and continues at a slowly- 
reducing rate to the present day. Our default models are 
truncation-only for simplicity; we show later that bursts of 
star formation on truncation affect the results very little. We 
assign an early-type galaxy a metallicity drawn from a dis- 
tribution with cT[ Fe / H ] =0.1 and a mean of sola r metallic itJEi 
motivated by the metallicity scatter inferred by Gallazz fet al.l 
(120061) . We use the multi-metallicity PEGASE stellar popu- 
lation model to explore the color evolution of the early-type 
galaxy population. We neglect, for the sake of simplicity, pos- 
sible differences in formation history as a function of galaxy 
masf3, while noting that the mean mass of our sample is 



8 We are assuming in the terminology of Schweizer & Seitzer (1992) that 
r 2 <lGyr. 

9 Galaxies with signatures of no current star formation but substantial star 
formation < 1 Gyr ago. 

10 such bimodality is washed out if transitions from blue to red take any 
longer than 1 Gyr 

1 1 The mean metallicity is of no importance in this paper, because we are 
exploring residuals from the color-magnitude relation. 

12 We see no significant difference in the behavior of this sample, when 
splitting into a high-mass and low-mass subsample, justifying this oversim- 



somewhat less than 1O U M . 

5.1.2. Major merger model 

The above model incorporates the effect of only one ma- 
jor merger: the one that created the early-type galaxy. The 
results of such a model are shown in Fig. [131 the results of 
a very similar model where 10% o f the stars are form ed in 
a burst on merging (motivated by iRobaina et al.l |2009|) are 
shown in Fig. [141 and are very similar to the truncation-only 
model. The upper left panel shows the build-up in the num- 
ber density of model red sequence, early-type galaxies, com- 
pared with observations of the build-up o f the red sequence 
stellar mass density from lBorch et al.l (|2006|) . The upper right 
panel shows the predicted evolution in the CMR zero point 
and scatter (think of this plot as a slice through the evolving 
color-magnitude relation at constant stellar mass): it is clear 
that there is a 'core' of red sequence galaxies, with a tail of 
blue early-type galaxies. 

In the center left panel, we explore the relationship be- 
tween truncation time (when star formation ceased) and color 
offset from the z = 0, t =12 Gyr red sequence. One can see 
that after truncation of star formation, the model galaxy red- 
dens rapidly, within ~ 1 Gyr, to within 0.1 mag of the red 
sequence. The subsequent reddening of the population is 
much slower. 'Typical' structured early-type galaxies, with 
A(g- r) ~ -0.04, would be interpreted to be between 2 and 8 
Gyr from their truncation event (depending on metallicity Q 
It is to be noted that this time scale (in its value and scat - 
ter) is in rough agreement with ISchweizer & S eitzer (1992); 
early-type galaxies less than 2a from the locus of the color- 
magnitude are expected to be more than 1 Gyr from their 'cre- 
ation event' . The oldest red sequence galaxies are redder than 
the core of the CMR. This is one of the reasons that the locus 
of the CMR appears to evolve slower than an ageing stellar 
population in such models; the average color of an early-type 
galaxy is always being pulled bluewards by new arrivals to the 
red sequence. 

We presented this relationship between truncation time 
and A(g - r), but can it have any relationship to real- 
ity (this, at some level, is the part of the question that 
ISchweizer & Seitzer] 119921 were less able to address quanti- 
tatively)? Put differently, in a model which is constrained to 
reproduce the number density evolution of early-type galax- 
ies, is the color distribution of present-day early-type galaxies 
correctly predicted? The result of this exercise is shown in 
the lower-left panel of Fig. \13\ The solid histogram shows 
the color offsets from the red sequence of the data, and the 
dashed histogram the major merger model. The two distribu- 
tions are clearly very similar: a KS test yields a > 10% chance 
of the two distributions being drawn from the same distribu- 
tion, which is a remarkable achievement for what is essen- 
tially a blind prediction, using only the observed red sequence 
number density growth and an assumed metallicity scatter. 

In the center and lower right-hand panels, we show light- 
weig hted ages of the sy stems (comparable to the ages derived 
by iGallazzi et all l2005) as a function of truncation time (cen- 
tral panel), and the histogram of light- weighted ages (lower 
panel). It is clear that the luminosity- weighted ages are 
skewed considerably towards older ages than the truncation 

plification. 

13 Recall that many of the early-type galaxies with slightly blue colors 
have little asymmetry in Fig. [6] and the large magnitude of this metallicity - 
dependent scatter in truncation time for a given color offset presents a natural 
and appealing explanation for this phenomenology. 
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Fig. 13. — A model in which major mergers create the early-type galaxy population. The top left panel s hows the creation rate of early-type (red, in this 
case) galaxies as a function of redshift from the models (solid line), and as measured by Borch et al. ( 2006). The top left panel shows the evolution of the 
predicted evolution of the color distribution of such a population as a function of the age of the population (the range 4- 12Gyr corresponds roughly to the 
interval 1 > z > 0). The middle left panel shows the relationship between color offset A(g-r) and the truncation time, and the middle right panel between 
luminosity-weighted age and truncation time. The bottom panels show the predicted (dashed lines) distributions of color offset and luminosity-weighted ages, 
and the observations (solid lines). It is clear that the present-day color distribution of early-type galaxies is reproduced by a model in which the early-type galaxy 
population is being built up at the observed rate, and that in such a picture there is a broad but scattered correlation between merger=truncation time and color 
offset from the color-magnitude relation. 



The merger origin of early-type galaxies 




Fig. 14. — A major merger model, but with a 10% burst of star formation just prior to the truncation of star formation. Our modeling conclusions are 
a burst of star formation on the final merger of a galaxy. 
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times, as is expected from any model galaxies that have a con- 
stant star formation rate before the truncation of star forma- 
tion. The distribution of ages is reasonably encouraging, and 
has the asymmetry characteristic of the data. The average age 
of the model is offset from the data; inasmuch as luminosity- 
weighted ages are considerably more model-dependent than 
color this offset could be relatively easily remedied (but is not 
here, for the sake of simplicity). For example, one could bring 
the start of star formation to ~ 10 Gyr in the past and adopt a 
different parameterization of the early stages of truncation of 
star formation in early-type galaxy progenitors. 

5.1.3. Minor merger model 

Of course, minor mergers are much more frequent than ma- 
jor mergers, and are increasingly thought to be an important 
feature of the evolution of the early-type galaxy population 
(see, e.g.. iNaab. Johansson & Ostriker 2009, Bezanson et al.l 
2009, Hopkins7taDl2010b|). The primary relevant effect here 
is that the minor merger of a gas-rich satellite onto a pre- 
existing early-type galaxy may drive the galaxy's color blue- 
wards for a short time and induce asymmetry. In Fig. [T5] 
we show a model that is identical to the major merger model 
above, but with the early-type galaxy population undergoing 
the accretion of a star-forming galaxy in the last 8 Gyr. The 
probability of accretion is modelled as being constant over the 
last 8 Gyr, with a total probability of 1 that a given galaxy ac- 
cretes a star-forming satellite over that period. The mass ratio 
is randomly chosen from a uniform distribution over the inter- 
val [0.0,0.4]; the metallicity of the satellite is assumed to be 
identical to the primary galaxy for modeling convenience. 

Fig. [15] shows many of the same phenomenologies as were 
seen in the major merger case (Fig. [13]). The evolution of 
number of red sequence galaxies is similar, as is the color 
distribution produced at to = 12 Gyr. The distribution of light- 
weighted ages is improved somewhat by the contribution of 
a small number of (bright) younger stars. Furthermore, the 
relationship between A(g - r) and the 'last' truncation time 
(the time of the last event to take place, be it a major or minor 
merger) is skewed considerably towards younger ages; it is 
already possible to be very close to the red sequence locus 
only ~ 1 Gyr after a minor merger. 

5.1.4. Bringing it together: our interpretation of the 
asymmetry-color relation 

Thus, the two key conclusions of the modeling effort are: i) 
the full color distribution of early-type galaxies at the present 
day is consistent with being built-up through truncation of 
star formation (as part of a merger event) at the rate inferred 
through observations of the build-up of the red sequence, and 
ii) minor mergers can lead to considerable short-lived color 
offsets (and presumably asymmetries), complicating consid- 
erably their interpretation. 

Despite the substantial scatter one expects in a major+minor 
merger model, the correlations shown in the central panels 
of Fig. [15] are still substantially stronger than the observed 
correlations be tween color/age and asymmetry (although the 
correlations of iTal et ail 1200 9 are almost as strong as those 
in the central panels of Fig. [15J it is possible that for ellip- 
tical galaxies residuals measured using deep data are well- 
correlated with time since last interaction). We quantitatively 
illustrate this issue in Fig. [TSJ We show the quantity ^ ast , the 
lookback time at which the last major or minor merger oc- 
curred in the major+minor mer ger model, statistical ly inferred 
from luminosity- weighted age (Galla zzi et al.l2005D . For each 



galaxy, a model galaxy with very similar luminosity-weighted 
age is found, and its ^ ast is assigned to the galaxy of interesQ 
Such an inferred ^ ast is purely statistical, by design reproduc- 
ing the trend and scatter in the middle right panel of Fig. IT31 
such a statistical estimate is for illustrative purposes only and 
does not represent an estimate of ^ ast that would be accurate 
in a galaxy-by-galaxy sense. We show the ^ ast inferred in this 
way as a function of asymmetry; the inset panels show the t\ ast 
distribution inferred for symmetric (A < 0.003) and structured 
(A > 0.003) galaxies. It is clear that asymmetry and ^ ast corre- 
late only very weakly; while structured galaxies have a weak 
tendency towards lower ^ ast , there are a number of structured 
galaxies with high t\ ast (i.e., rather red stellar populations). Put 
differently, this figure attempts to encapsulate quantitatively 
the fact that neither the stellar population parameters nor the 
asymmetry are perfect clocks: the scatter in star formation 
history, gas content, mass ratio and details of the merger or- 
bits appear to add very significant scatter to the relationship 
between the stellar population and structural 'clocks'. 

In this kind of picture, while one expects a broad correla- 
tion between color offset from the red sequence and asym- 
metry (as is observed), there will be a lot of scatter in this 
relationship for astrophysical reasons. Some relatively re- 
cent minor accretions, or interactions bet ween already non- 
star-forming progenitors (dry merge rs; Ivan Dokkuml [2005: 
iBell et al.l 2006: M cintosh et ai1l2008l) . might be richly struc- 
tured but already red; whereas some more ancient major in- 
teractions may remain blue for much longer, showing little 
asymmetry but blue colors. In such a picture, finding galaxies 
to fill the 'King gap' is not a particularly well-posed exer- 
cise. It is of course possible to define samples that will likely 
contain some fraction of the most recently-formed early-type 
galaxies. Such a candidate sample is presented in Fig. [17] and 
Table [2] where we present the 20 bluest objects in our sam- 
ple, with A(g-r) < -0.14. Inspection of Fig. [T71 shows that 
many of them clearly show asymmetries generated by tidal ef- 
fects; yet, some of this sample are not obviously asymmetric, 
illustrating the challenges of uniquely identifying a sample of 
'young' early-type galaxies. 

Yet, we would argue that there should be little concern 
about this difficulty to uniquely identify intermediate- age 
merger remnants. Much more important in our view is the 
finding that in this model, in which early-type galaxies are 
being produced at the right rate by major interactions (with or 
without the additional accretion of satellites), the correct dis- 
tribution of early-type galaxy colors is straightforwardly re- 
produced. This indicates that as an ensemble, the (stellar pop- 
ulation) properties of the early-type galaxy population are in 
accord with a model in which they are being built up by merg- 
ing at the observed rate. Our result confirms a broad link be- 
tween asymmetric (largely tidally-induced) structure and stel- 
lar population parameters, confirming a qualitative prediction 
of the merger hypothesis. Detailed Af-body and hydrodynam- 
ical modeling will be necessary to test this link quantitatively, 
by providing directly the predicted distribution of galaxies in 
the color/age vs. asymmetry plane as predicted by assembly 
histories derived from simulations of galaxy formation in a 
cosmological context. 



14 By doing this, we end up with more scatter than is necessary in the final 
fiast -asymmetry relation. We choose a model galaxy with a given luminosity- 
weighted age at random, whereas nature is likely to ensure that highly struc- 
tured galaxies tend to have a lower than average t\ ast . 
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Fig. 15. — A model in which major mergers create the early-type galaxy population, with the addition of minor mergers. The figure is formatted similarly to 
Fig- El except that in this case the middle panels show population parameters as a function of the truncation time of the last event to take place, be it a major 
or a minor merger. Again, it is clear that the color distribution of early-type galaxies is well-reproduced by such a model, but that the addition of minor mergers 
scatters further the already scattered correlation between color offset from the CMR and last truncation time. 



5.2. Musings on the role of dissipationless (dry) merging 

In recent years, it has been demonstrated that largely dissi- 
pationless merging between already-formed early-type galax- 
ies (dry mergers) play a signi ficant role in the build-up of the 
early- t ype galaxy popu lation (Bell et al. 2004a; van Dokkum 
[20051: iBeUeTa l. 2006; lLotz et al.l 120081: iFaber et alT l2007: 
ISkelton et al.l 120091) . In particular, the most massive 
early-type galaxies appear to grow only by dry me rging 
(IBell et al.ll2004aHFaber et al.l2007tlMcTntosh et al.ll2008h . al- 
though what that growth rate is remains somewhat unclear 



(Ivan Dokkuml l2005[ IBell et al.l 120061: iMasiedi et al.l 2006; 
Scarl ata et alll2Q07b iMdntosh et al.ll2008l) . 

In this context, it is of interest to note the population of 
relatively red ( A(g - r) ~ 0) and richly structured E and SO 
galaxies in Figs. [5] and [TT] We show a set of 20 randomly 
selected red highly asymmetric galaxies with A(g-r) > -0.1 
and A > 0.02 in Fig. Q21 These systems are clearly undergo- 
ing (or have undergone) interactions, both minor and major, 
with already early-type galaxies. The relative prominence of 
this tail of red asymmetric E/S0 galaxies appears to be rela- 



16 



Gyory & Bell 



TABLE 2 

Parameters of the 20 bluest E/SO galaxies 



RA 


Dec 


A(g-r) 


log (age [Gyr]) 


log (Z/Z ) 


A 


224.1003 


49.6960 


-0.14 






0.002 


166.3744 


3.6600 


-0.14 


9.56 


-0.28 


0.003 


142.5555 


49.4882 


-0.14 






0.005 


227.0540 


56.4228 


-0.14 


9.32 


-0.03 


0.012 


169.7402 


58.0566 


-0.15 






0.003 


17.9440 


-0.6645 


-0.15 






0.002 


226.6311 


40.6962 


-0.16 






0.013 


202.2038 


53.4430 


-0.18 






0.009 


135.1501 


46.6863 


-0.19 






0.007 


167.6620 


3.7555 


-0.19 


9.51 


-0.46 


0.005 


32.8899 


13.9171 


-0.19 


8.84 


-0.04 


0.016 


11.2468 


-8.8897 


-0.20 


9.49 


0.03 


0.024 


50.8886 


-0.4385 


-0.20 






0.002 


135.7590 


40.4340 


-0.21 






0.012 


153.9239 


7.0522 


-0.21 






0.013 


199.9037 


3.0327 


-0.22 


8.87 


-0.20 


0.037 


129.9133 


3.8285 


-0.24 






0.022 


210.1692 


-1.9217 


-0.29 






0.004 


120.8668 


25.1026 


-0.31 






0.039 


205.3000 


1.7798 


-0.42 


9.39 


-1.08 


0.013 



to properly interpret this behavior. Two aspects will clearly 
play a role in setting the relative prominence of blue vs. red 
richly- structured early-type galaxies: the ratio of blue vs. red 
merger partners typical of the mergers that create early-type 
galaxies, and the issue of when one classifies a merger rem- 
nant as an early-type galaxy (almost always for two interact- 
ing early-type galaxies, as opposed to only at very late times 
for the remnant of a merger between two gas-rich galaxies). 

It is also of historical interest to ask why there were rela- 
tively few such systems in the work of Schweizer & Seitzer 
(1992). Number statistics clearly will play a role, but also im- 
portant is the differing natur e of our measures of t idal struc- 
ture in early-type galaxies. ISchweizer & Seitzel (Il992|) at- 
tempted to measure fine structure characteristic of mergers 
between gas-rich systems, whereas we explore asymmetry, 
which is also sensitive to large-scale asymmetries and de- 
bris fields, more characteristic of mergers between early-type 
galaxies. This aspect, along with details of initial sample se- 
lection, is likely to play a role in this small difference between 
our two studies. 

5.3. Comparison with Tal et ah (2009) 

iTal et alJ (120091) choose a complementary method to ex- 
plore this issue, by quantifying the residuals from a smooth 
model fit for a sample of elliptical galaxies and correlat- 
ing them with color offset from the CMR (residuals from a 
smooth fit for SOs would reveal the disk/bulge transition and 
bars, and would be significantly less meaningful). They used 
very deep data for 55 systems and traced tidal structure out 
to large radius, where dynamical times are long. They found 
that a large fraction of systems were disturbed, and found a 
clear but scattered correlation betwe en color offset f r om th e 
CMR and residuals, as seen by Sch weizer & S eitzer (1992). 
Their results can be regarded as a high sensitivity analysis 
of the main body of the correlation shown in Fig. [6l where 
their higher sensitivity to faint debris has brought the color- 
residual correlation into sharp focus. Our result complements 
this analysis by i) using a larger sample and including SO 
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Fig. 16. — Time of last major or mino r merger, i nferred statistically from 
the luminosity-weighted age (from Gallazzi et al. 2005), as a function of 
asymmetry. The dashed line shows the dividing line between 'symmetric' 
galaxies with A < 0.003 and 'structured' galaxies with A > 0.003. The inset 
panel shows the histograms of inferred last merger times for the symmetric 
and structured samples. 

tively independent of luminosit}0 (somewhat in contrast to 
the naive expectation that dry interactions/mergers are more 
important for more luminous systems), indicating that such 
interactions are possible for systems of a wide range of lu- 
minosities from L* /5 right up to the most luminous systems 
probed in our study ~2- 3L*. 

It is of interest that many of the most structured early-type 
galaxies have this signature. Without a wide range of TV- 
body /hydrodynamical simulations to guide us, it is difficult 

15 This was checked by splitting the sample in half by luminosity. There 
was no significant difference between the red asymmetric early-type galaxy 
fractions of the bright and faint bins. 
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Fig. 17. — The 20 bluest E/SO galaxies. Each box is approximately 55 x 55 arcsec. A lOkpc horizontal line is shown for each galaxy. 



galaxies, showing it to be generally applicable to the popu- 
lation as a whole (albeit with a lower sensitivity measure of 
structure, largely because of limitations in SDSS depth), and 
ii) using age measures to explicitly illustrate that the main 
driver of such behavior is an age-structure relation. 

5.4. Limitations of this approach, and outlook 

As alluded to in the previous sections, a significant (and 
largely unavoidable) limitation of our study is the choice of 
asymmetry as a metric of tidally-induced disturbances. Our 
choice of asymmetry was motivated by the overriding consid- 
eration that the measurement be algorithmic and reproducible; 
asymmetry was chosen over other reproducible descriptions 
of structure (such as resi duals from a sm ooth model fit; e.g., 
iMcIntosh etaLll2004l and lTal et aLll2009b to avoid symmetric 
bars or rings from contributing to our measure of structure. 
Yet, a key limitation of our metric is that systems can be asym- 
metric for a variety of reasons: asymmetry may be produced 
by fly-bys or mergers, major or minor interactions, ongoing 
or past interactions, gas-rich or gas-poor mergers. As such, 
asymmetry is a valuable but ultimately blunt tool; in this re- 
spect, our work is a useful way to start attacking the problem, 
but is not capable of answering higher-level, more detailed 
questions about the origin of the early-type galaxy population. 

In principle, one could attempt to construct a richer descrip- 
tion of structure, which would be more clearly correlated with 
particular types of interaction (e.g., a measure of the spatial 
scale of tidal debris, or a measure of the luminosity in tidal 
debris). Calibration with simulations would be a critical step 



in the establishment of such metrics (indeed, we have not yet 
calibrated our measure of asymmetry as part of this work). 
Yet, in practice, such metrics have prov en difficult to reliably 
constr uct; the fine structure parameter of lSchweizer & Seitzeri 
(1992) is a good example of such a metric, which may be more 
physically meaningful but is difficult to reproduce. 

6. CONCLUSIONS 

One of the key predictions of the merger hypothesis for the 
origin of early-type galaxies is that tidally-induced asymmet- 
ric structure should correlate, at least in a broad sense, with 
signatures of a rela tively young stellar popula tion. Motivated 
by this argument, ISchweizer & SeitzeT(ll992|) studied the re- 
lationship between color offset from the CMR and tidally- 
induced fine structure, finding a correlation between the two 
quantities. They modeled the colors of these structured early- 
type galaxies, rinding that such colors can be reproduced by 
a variety of models in which star formation truncates reason- 
ably quickly, with typical timescales of a few Gyr to reach the 
colors characteristic of ancient early-type galaxies (although 
with much model dependence). 

In this paper, we have re-examined this issue, incor- 
porating a number of im provements over the work of 
ISchweizer & Seitzed (119921) : a sample from 2M ASS/NED 
and the SDSS that is nearly ten times larger with well- 
measured colors; the use of a non-ideal, but reproducible and 
meaningful metric for tidally-induced structure, asymmetry; 
the use of luminosity-weight ed ages and metallicities from the 
work of Galla zzi et al.l {2005) using the SDSS; and, the ability 
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Fig. 18. — A sample of 20 red, asymmetric E/SO galaxies. Each box is approximately 55 x 55 arcsec. A lOkpc horizontal line shown for each galaxy. 



to construct a more physically-motivated model of early-type 
galaxy evolution with which to i nterpret the results. 

We found, i n agre ement with ISchweizer & Seitzeri (|1992|) 
and iTal et al.1 (|2009|) . a correlation between offset from the 
CMR and asymmetry. Inspection of asymmetric systems gave 
weight to the notion that the main driver of asymmetries was 
galaxy interactions (although a variety of interactions can cre- 
ate asymmetries, leading to a natural ambiguity when inter- 
preting the results). We demonstrated, for the first time, that 
age effects are driving this correlation. A fraction of asym- 
metric early-type galaxies have normal colors/ages, character- 
istic of mergers between already-formed early-type galaxies 
(dry mergers). The empirical correlation between stellar pop- 
ulation age and tidally-induced asymmetries, for the bulk of 
the population, is the key result of this paper, and is consistent 
with the basic prediction of the origin of early-type galaxies 
through galaxy merging, and their modification by accretion 
of gas-rich satellites. 

We interpreted these results in the context of a model in 
which the number of early-type galaxies is increasing con- 
stantly with time, in quantitative agreement with measure- 
ments of the evolving number density of red sequence (early- 
type) galaxies. Assuming that the ultimate effect of a galaxy 
merger is to truncate star formation (through either gas con- 
sumption or AGN feedback), we find that such a model repro- 
duces the distribution of color offsets from the CMR; put dif- 
ferently, the merger hypothesis appears to give the correct dis- 
tribution of present-day early-type galaxy colors. Under the 



assumption that both light- weighted stellar age and asymme- 
try are good 'clocks' of the time elapsed since the last major or 
minor merger, such a model predicts a rather tighter correla- 
tion between age and asymmetry than is observed. This indi- 
cates that astrophysical sources of scatter, e.g., the star forma- 
tion histories of the progenitors, their gas content, and mass 
ratios, play an important role in setting the properties of the 
remnants. Thus, the broad correlation between age and asym- 
metry is in qualitative (but not yet quantitative) agreement 
with such a model. The existence of such a large scatter pre- 
cludes the unambiguous identification of 'young' early-type 
galaxies from stellar population and asymmetry data alone; 
confirming the merger hypothesis using a few individual ex- 
amples will be a challenging and ambiguous exercise. 

We conclude that the properties of the early-type galaxy 
population as a whole are consistent with the bulk of them 
being formed by galaxy merging (some of these mergers are 
between already-formed early-type galaxies). In the case of 
stellar population parameters, there is a quantitative match 
between toy models of merger-driven growth and the obser- 
vations. In the case of asymmetries, the existence of a broad 
and scattered correlation between stellar population parame- 
ters and asymmetry is in qualitative agreement with the expec- 
tations from the merger hypothesis, but our toy models (and 
even more sophisticated models) are currently incapable of 
providing expectations for how asymmetry should vary with 
time and interaction phase. Further work, using simulations 
to predict and calibrate the distributions of asymmetries, and 
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using more sophisticated descriptors of galaxy structure, will 
help to more deeply explore this issue in the next years. 
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